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ABSTRACT 

We report the detection of radio emission and orbital motion from the nearby star-brown dwarf binary 
WISE J072003.20—084651.2AB. Radio observations across the 4.5-6.5 GHz band with the Very Large 
Array identify at the position of the system quiescent emission with a flux density of 15±3 /xJy, and 
a highly-polarized radio source that underwent a 2-3 min burst with peak flux density 300±90 /iJy. 

The latter emission is likely a low-level magnetic flare similar to optical flares previously observed 
for this source. No outbursts were detected in separate narrow-band Ha monitoring observations. 

We report new high-resolution imaging and spectroscopic observations that confirm the presence of a 
co-moving T5.5 secondary and provide the first indications of three-dimensional orbital motion. We 
used these data to revise our estimates for the orbital period (d.lfl^'g yr) and tightly constrain the 

orbital inclination to be nearly edge-on (93?6'^^i^), although robust measures of the component and 
system masses will require further monitoring. The inferred orbital motion does not change the high 
likelihood that this radio-emitting very low-mass binary made a close pass to the Sun in the past 
100 kyr. 

Keywords: binaries: visual — stars: chromospheres — stars: individual (WISE J072003.20—084651.2) 

— stars: low mass, brown dwarfs — stars: magnetic field 


1. INTRODUCTION 

WISE J072003.20-084651.2 (hereafter 

WISE J0720 — 0846) is an M9.5 dwarf originally 
identified by iScholzl (|2014ll as a possible new mem¬ 
ber of the local 8 pc sample. It was previously 
missed in searches for nearby low-mass dwarfs due 
to its low Galactic latitude {b = 2.3°) and small 
proper motion (121.7±0.3 mas yr~^). Su bs equen t 
follow-up observation s by IKirkoatrick et al.l (1201411: 
iBurgasser et al.l (|2015L hereafter BIS'): and llvanov et al.l 
(|2015L hereafter 115) confirmed the late-type nature 
and proximity of the source (6.0±1.0 pc), and have 
identified it as having an unusually high recessional 
velocity (-1-83.8 k m s~^), indicat i ve of old disk/thick 
disk kinematics. IMamaiek et al.l (|2015ll have deduced 
that WISE J0720—0846 had one of the closest stellar 
approaches to the Sun inferred to date, passing within 
0.25^0.07 pc over a period of 60,000-85,000 years ago, 
possibly penetrating the outer Oort Cloud. 

Evidence for a T-type brown dwarf companion to 
WISE J0720—0846 was reported by B15, based on both 
near-infrared spectral analysis and the presence of a 
candidate resolved source in high-resolution laser guide 
star adaptive optics (LGSAO) imaging. In support of 
this, 115 reported infrared excess in 11 and 22 fim pho- 
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tometry from the Wide-he ld Infrared Survey Explorer 
fWISE: IWright et al.ll201fll l that could be attributed to 
a low-temperature secondary. However, the detection of 
the secondary was marginal due to its close separation 
(140 mas), large magnitude difference with the primary 
{AH = 4), and poor LGSAO correction, and neither B15 
nor II5 detected significant radial velocity (RV) variabil¬ 
ity over overlapping ^3 month timescales. Confirming 
the presence of this putative companion and assessing 
the degree of its gravitational perturbation on the pri¬ 
mary is important for determining an accurate parallax 
and space motion for WISE J0720—0846. This system 
is only the second (candidate) binary to be identified 
among the 14 known late-M dwarfs within 10 pc of the 
Sun. 

Another remarkable trait of WISE J0720—0846 re¬ 
ported in B15 and 115 is its weak yet highly variable 
magnetic emission. Quiescent Ha emission was observed 
to be at or below typical values for M9-L0 dwarfs, with 
a relative power of log^o Lna/Lboi ~ —5. However, B15 
reported the detection of multiple flaring events, both 
in white-light photometry and Ha line emission. These 
flares were infrequent (<1% effective duty cycle) but 
nevertheless produced order-of-magnitude variations in 
nonthermal emission. Unfortunately, the X-ray limit of 
log;^g Lx/Lboi < —3.2 reported by 115 does not provide a 
stringent constraint on the high-energy nonthermal emis¬ 
sion from this source. Overall, WISE J0720—0846 ap¬ 
pears to be similar in magnetic behavior to the rapidly 
rotating BRI 0021-0040 (iBasri fc Marcvll995HReid et al.l 
119^ . being weakly active but with occasional strong 
bursts. Deducing how these bursts relate to age or rota¬ 
tion, or in the case of WISE J0720—0846 interaction with 
a (putative) close companion, may provide critical clues 
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for understanding the overall decline in optical and X-ray 
magnetic ac t ivity a cross the M dwarf/L dwarf transition 


(IGizis et al.l 120001: 

Kirknatrick et al.l 120001: IWest et al.l 

12004 IStelzer et al.l 

20061: iSchmidt et al.l 120141 and con- 

trarv trends in nonthermal radio emission (Bere'eill2006l: 


iMcLean et alJ[2?)T^ lAntonova et aljr2013ll . 


In this article, we report new observations of 
WISE J0720—0846 at optical, near-infrared and radio 
wavelengths which confirm both the bursting and binary 
nature of this source. In Section 2 we report the de¬ 
tection of low-level quiescent and bursting radio emis¬ 
sion based on data obtained with the Karl G. Jansky 
Very Large Array (hereafter VLA), and limits on Ha 
variability from asynchronous narrow-band photometric 
monitoring. In Section 3 we report new imaging and 
high resolution spectroscopic observations that confirm 
the presence of a T5.5 companion, and provide the first 
indications of orbital motion. We use these data to make 
constraints on the orbital configuration of the system in 
Section 4. Our results are discussed in Section 5. 

2. MAGNETIC EMISSION FROM WISE J0720-0846 

2.1. Observations 

WISE J0720—0846 was observed with the VLA 
(project code I4B—313, PI Burgasser) in the com¬ 
pact C-configuration (baselines of 0.035—3.4km) on 2014 
November 11 from UT 08:48:50 to 14:47:48. The WIDAR 
correlator was set up for C-band continuum observa¬ 
tions with two basebands, each having eight 128 MHz 
sub-bands centered at 5.0 GHz and 6.0 GHz, for a fre¬ 
quency range of 4.488 GHz to 6.512 GHz and a total 
bandwidth of ~2 GHz. Each sub-band had four polar¬ 
ization products (RR, LL, RL, LR) and sixty-four 2 MHz 
channels; a 5 sec dump rate was used. After observa¬ 
tion of our primary calibration source 3C48 to set the 
absolute flux scale and measure the complex bandpass, 
we conducted a sequence of 4 min cycles, with 3 min 
on WISE J0720—0846 and 1 min on the gain calibrator 
QSO J0730—1141. This observational strategy optimizes 
image quality by frequent monitoring (and hence correc¬ 
tion) of phase fluctuations. 

All data were reduced with the A stronomical Im age 
Processing Software package (AIPS; iGreisenI 1200311 fol¬ 
lowing best practices for wide-band radio data reduction. 
Radio-frequency interference (RFl) was present through¬ 
out the observation, being especially persistent and 
strong in the sub-bands covering 5.782 GHz to 6.512 GHz, 
and also present at varying levels in several other sub¬ 
bands. All data affected by RFl were removed. In ad¬ 
dition, sub-band 0 of each baseband failed to produce 
robust data, and these measurements were also removed. 
After data flagging, our average frequency was 5.27 GHz 
and our total bandwidth was reduced to si 1.2 GHz. 

2.2. Quiescent Emission 

The final cleaned and calibrated data were analyzed 
by imaging and performing time-series analysis on the 
uu-data. Figure [T] shows imaging of the entire data 
set (all unflagged times and frequencies) at the coordi¬ 
nates of WISE J0720—0846, which reveals a weak, elon¬ 
gated source. This is interpretted to be a blend be¬ 
tween weak quiescent emission from WISE J0720—0846 
and an unrelated faint background source, each having 
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Figure 1. VLA image of the WISE J0720—0846 field integrated 
over 4.5—6.5 GHz and over the entire time series on 2014 Novem¬ 
ber 14 (UT). Contour levels in flux density are labeled at bot¬ 
tom. The beam shape is indicated in the lower left corner and 
has dimensions of 478x375 with position angle of 11?5. Emission 
from WISE J0720—0846 emerges from an extended source span¬ 
ning ~10'' along a north-south axis, likely arising from combined 
emission from the target and an unassociated background source. 

similar flux densities of ~15/xJy0 Using the AIPS task 
DFTP L and the procedures described in lOsten fc WolkI 
(Hooi), we analyzed the time-series data at the peak 
flux position, shown in Figure [21 This series reveals 
weak but nonzero emission over the entire course of 
the observation, with a short burst of emission at UT 
13:21:37. Masking out the burst, we measure a mean 
quiescent flux density of 15±3 ^Jy with no evidence of 
statistically significant variabilityQ This is the weakest 
quiescent emission detected for a very low mass dwarf 
to date, comparable to the quasi-quiescent emission re- 
ported for the T6.5 d warf 2MASS J10475385-b2124234 
(|Williams et al.l I2013L 16±5 /iJy). The apparent flux 
density translates into a specific radio luminosity of 
\ogiQLi, = (7±2)xl0^^ erg s“^ Hz“^ at the 6.0±1.0 pc 
distance of WISE J0720—0846, and log^g Vp^L^/Lboi = 
—8.49±0.19 assuming log^g Lboi/LQ = —3.60±0.05 (B15; 
115) and Vpk = 5.28 GHz. The specific luminosity is 1- 
2 orders of magnitude below measurements of previo usly 
detected late M and L dwarfs ([Antonova et al.l[2013f) and 
the relati ve luminosity falls be low the activity-rotation 
relation of lMcLean et al.l (120121 1 given this source’s v sin i 
= 8.6±0.8 km s“^ (Section 3.4). Indeed, in terms of rota¬ 
tion and radio p ower, WISE J0720—0846 is s imilar to the 
M7 dwarf VB 8 ([Krishna murthi et al.||19 99ll and the M 8 
dwarf DENIS J1048-3956 (iBurgasser fc PutmanI [200511 . 
both nearby late M dwarfs with weak but variable mag¬ 
netic emission and modest rotation rates. The weakness 
of the quiescent emission prevents us from determining 
either its spectral behavior across the 4.4-6.5 GHz band 
or its polarization. 

2.3. Bursting Emission 

® A resolved-object fit to this elongated source yields a total 
integrated flux density of 28±8 /rjy; the map rms noise level is 
2.8 /rjy bm“^ 

® For 1 minute sampling, we measure = 296.0 for 295 degrees 
of freedom, yielding a p-value of 0.46. Note that with the bursting 
emission included = 382.5 and p <0.1%. 
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Figure 2. (Top) Time series of radio fiux from WISE J0720—0846 in the 4.5-6.5 GHz band with 1 min sampling; uncertainties are 
indicated by error bars. The red dashed and dotted lines indicate the mean flux level and Icr uncertainty, after excluding the bursting 
emission at UT 13:21. (Bottom) Close-up view of the burst emission, with data sampled every 10 s. Breaks correspond to observations of 
the secondary calibrator. The best-fit Gaussian (red solid curve) and exponential decay (blue curves: dashed for peak emission at the start 
of calibration observation, solid for peak emission at the end) models are overplotted. 
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Figure 3. VLA image of the WISE J0720—0846 field in Stokes I (left) and V (right) integrated over 4.5-6.5 GHz and over the period 
UT 13:19:07 to 13:24:07. The burst emission seen in the time series data is associated with a point source at the expected position of 
WISE J0720—0846. The Stokes V image shows a negative source of comparable brightness indicating nearly 100% left circular polarization. 


A close-up view of the radio burst around UT 13:21 is 
shown in Figure [5] The burst appears to have begun dur¬ 
ing observations of the secondary calibrator, so our data 
sample it only at or after its peak emission. Reimaging 
of the radio data in the period around this burst shows 
it to be associated with a bright point source at the co¬ 
ordinates expected for the brown dwarf (Figure [3]) Q The 
burst is highly polarized, with emission in Stokes / of 
175±28 /iJy and Stokes V of —150±28 fiJy, implying 
left-handed circular polarization of 84^20% (1 ^ equiva¬ 
lent uncertainties); i.e., consistent with full polarization. 
There is no evidence of spectral variation in the burst to 
within our measurement uncertainties. 

We modeled the time series emission as both a ro¬ 
tating spot (“pulse”) feature with a Gaussian profile, 
and as an exponentially declining flare, oc e“*/'^ af¬ 
ter peak emission. The former yields a peak flux density 
of 310±40 ^Jy, ~20 times greater than the quiescent 
emission, and consistent with a total emitted energy of 

^ The source identified during the flare is free of contamination 
because the instantaneous sensitivity over the small time range 
imaged is insufficient to detect the background blend source. 


(1.3±0.4) X10^^ erg over the 4.5-6.5 GHz band (assum¬ 
ing a flat spectrum). The flare model yields a decay time 
constant of t = 3.5±0.9 min and, depending on whether 
the burst initiated at the start or end of the calibration 
period, a total emitted energy of (1-3) x 10^"^ erg. 

The burst emission occurs for at most 3 min during 
our 5 hr of observation of WISE J0720—0846, or 1% of 
the total time on-source. This is consistent with the flar¬ 
ing duty cycle inferred from the aperiodic white light 
bursts reported in B15. Thus, we favor an infrequent 
flare mechanism for this emission as well. We neverthe¬ 
less note that the pulse/spot model implies an emission 
region Ril-2° in longitudinal extent (assuming a rota¬ 
tion period of ~14 hr; see Section 4), whi ch is consistent 
with previous periodic pulse d etections (iHallinan et al.l 
120071 : TWilliams fc Bergeil l2015f) . Longer-term monitor¬ 
ing would be required to distinguish between the flare 
and pulse hypotheses. 

2.4. Ha Monitoring 

WISE J0720—0846 was monitored over four nights on 
2014 February 18, 20, 22 and 24 (UT) using the facility 
CCD camera on the 1-meter Nickel telescope at Lick Ob- 
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Table 1 

Ha Monitoring Observations 


UT Date 

MJD 

^int 

(s) 

Nobs 

UT Start 

Duration 

(hr) 

Mean S/N per 
Measurement 

2014 Feb 18 

56706.26 

900 

11 

4:42:32 

3.03 

4 

2014 Feb 20 

56708.24 

900 

10 

4:38:51 

2.63 

3 

2014 Feb 22 

56710.24 

300 

38 

3:55:13 

3.97 

4 

2014 Feb 24 

56712.28 

300 

13 

4:55:02 

3.52 

3 


servatory (Table [T|). The CCD was configured for 2x2 
binning for a pixel scale of 0"37 pixel”After acquisition 
and centering at /-band, WISE J0720—0846 was moni¬ 
tored without dithering through the narrow-band Ha fil¬ 
ter (Ac = 6557 A, AA = 15 A), with integration times 
of 900 s on February 18 and 20 and 300 s on February 
22 and 24. Total monitoring periods per night spanned 
2.63-3.97 hr, but due to overheads and pauses during oc¬ 
casional clouds the on-source time totalled 7.25 hr over 
the entire observing run. Bias frames and quartz flat 
field lamps were also acquired each night for detector 
calibration. 

Data were reduced using standard image reduction 
techniques. Aperture photometry of WISE J0720—0846 
and nearby non-saturated stars was measured using a 
variable aperture scaled to encapsulate 80% of each 
source’s peak brightness, and an annulus of 63-100 pixels 
(23"-36") was used to subtract foreground emission. We 
did not observe a photometric calibrator during these ob¬ 
servations, so no attempt was made to measure absolute 
Ha fluxes. Instead, we used the mean flux of non-variable 
stars in the field of view to compute a reference light 
curve, and used this to normalize the photometry for 
WISE J0720—0846 over the course of each night. Uncer¬ 
tainties were dominated by source photometry, of order 
20-50%. 

We found no significant flux variations over the course 
of the four nights of observation. At the measured noise 
level, we can rule out bursts 2-3 times above quiescent 
emission. B15 reported a single Ha flare with a line flux 
5 times greater than quiescent emission, and an order of 
magnitude brighter than the local continuum. Such a 
flare would have been easily detected in our observations 
had it occured. Assuming a flare period of ~5 min, we 
infer a flare indicidence rate of <1%, similar to the white 
light flare rates reported in BI5 and the radio flare re¬ 
ported here. All of these observations reinforce the con¬ 
clusion that WISE J0720—0846 is a weakly active and 
infrequently bursting source. 

3. CONFIRMING THE MULTIPLICITY OF WISE J0720-0846 

3.1. High Resolution Imaging 

WISE J0720—0846 was re-observed with the Near- 

InfraRed Came ra _ 2 (NIRC2) and LGSAO system 

(jvan Dam et al.l 120061 : IWizinowich et al.l 1200611 on the 
Keck H 10 m Telescope on 2015 January II (UT) in 
mostly clear and windy conditions with 078 seeing. The 
narrow field-of-view (FOV) camera was used to obtain 
dithered observations in broad-band MKC0 /, H and 
Ks filters, and the medium-band CH^s filter sampling 

^ Mauna Kea Observatories near-infra red filter set 
IlSimons fc Tokunagal 1^021 : ITokunaga et al]l2002l l. 


Table 2 

Resolved Photometry and Astrometry of WISE J0720—0846 


Parameter 

Value 

AJ 

2.92±0.07 

ACH4S 

3.05±0.08 

AH 

3.85±0.11 

AKs 

4.07±0.20 

p (mas) 

197±3 

P (AU) 

1.19±0.21 

PA (°) 

256.7±0.6 

Primary SpT'^ 

M9.5±0.5 

Secondary SpT'^ 

T5.5±0.5 


^ Based on spectral template fitting with templates con¬ 
strained to have the same relative flux scaling as measured in 
the JHKs NIRC2 bands. 


1.54-1.65 /tm. The R = 16.8 mag field star USNO 0812- 
0137390 was used to correct for tip-tilt aberrations. We 
achieved better Strehl ratios than observations reported 
in B15, 10-20% depending on wavelength. This allowed 
us to easily resolve the candidate companion reported in 
that study in all four filters, as shown in Figure SI 

3.2. Common Proper Motion and Orbital Motion 

To extract relative photometry and astrometry, we per¬ 
formed point source function (PSF) fitting of each indi¬ 
vidual image using a two-dimensional asymmetric Mof¬ 
fat profile optimized to the PSF of the primary com¬ 
ponent (i.e., with the secondary masked). Measure¬ 
ments are reported in Tabled] We find a separation of 
197±3 mas at position angle 256?7±0?6, wider than, and 
at a marginally distinct position angle as, the candidate 
source previously reported (139±14 mas at 262°±2°). As 
the source detected in these data would have been eas¬ 
ily resolved in prior observations, we conclude that sta¬ 
tistically significant relative motion has been observed 
between the two epochs. 

Relative motion can be due to differential motion be¬ 
tween two physically unrelated sources or orbital motion 
in a gravitationally bound binary. Figure [5] displays the 
estimated center of mass motion and component posi¬ 
tions between our 2014 and 2015 imaging data, assuming 
systemic astrometry from B15 and q — 0.4 (see below). 
It is clear that the predominant motion of both sources 
is co-aligned with the proper motion of the system, par¬ 
ticularly in declination; the change in position angle is 
inconsistent with the secondary being an unmoving back¬ 
ground source at the I6cr confidence level. We therefore 
determine that WISE J0720—0846 is a co-moving binary 
system, and identify significant astrometric orbital mo¬ 
tion over a one year-period. 

3.3. Improved Determination of the Secondary 
Classification 
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Figure 4. Resolved imaging of the WISE J0720—0846 system with NIRC2/LGSAO. Images are aligned with North up and East to the 
left, and each box displays an angular scale of 1". Images are logarithmically flux-scaled to make the faint secondary visible. 



Aa cos 5 (arcsecond) 

Figure 5. Positions of WISE J0720—0846A (red points) and B 
(blue points) in our 2014 and 2015 LGSAO images relative to 
center-of-mass astrometric motion (black points and dotted line) 
based on values from B15 and assuming q = 0.4. The 3cr uncer¬ 
tainties on the center-of-mass position at the 2015 epoch relative 
to 2014 are indicated: (Tc« = 4 mas, ag = 0.8 mas. The relative 
motion of the two sources is consistent with physical association 
and orbital motion. 

Our new observations provide greatly improved rel¬ 
ative photometry of the two components in the four 
bands measured, and refine the coarse AH estimate 


from B15 (Table [2]). There is a significant difference 
in relative brightnesses in the CH^s and 7J-band fil¬ 
ters consistent with strong CH 4 absorption in the sec¬ 
ondary. We used the relative photometry to better con¬ 
strain the spectral types of the components through spec- 
tr al template fittin g Foll owing the procedures described 
in iBurgasser et all (I201H ). we combined 512 M7-L1 and 
125 T0-T7 spectral templat es from the SpeX Prism Li¬ 
brary (SPL: lB7irgassei]l20l4) . scaled so that that relative 
spectrophotometry aneed to within Icr of the NIRC2 
JHKg measurements O We compared the 16394 binary 
templates that satisfied these constraints to the com- 
bined light SpeX spectra of WISE J0720—0846 from both 
iKirkoatrick et all (120141 ) and B15, following the methods 
described in the latter paper. Figure | 6 ] shows the best fit¬ 
ting template to the B15 spectrum. Both analyses yield 
identical results, with component types M9.5±0.5 and 
T5.5±0.5. The secondary is a half subtype later than, 
but formally consistent with, the classification reported 
in B15. 

3.4. High Resolution Spectroscopy 

New high resolution optical and near-infrared spec¬ 
troscopy of WISE J0720—0846 were obt ained with the 
Near InfraRed Spectrometer (NIRSPEC; iMcLean et al.l 
120001) on the Keck II telescope on 2014 Decem ber 8 (UTb 
and with the Hamilton echelle spectrograph (lVogtlll9^ 
on the Lick Observatory Shane 3 m telescope on 2015 
March 10 (UT). Data were acquired and reduced as de- 

® NIRC2 CH^s photometry was not used for this analysis as the 
filter profile was unavailable in a digital format (H. Tran, 2015, 
priv. comm.). 






























6 



Figure 6. Best-fit spectral binary template (purple line) to the 
combined-light SpeX spectrum of WISE J0720—0846 (black line) 
based on scaling spectral templates to the relative JHKs photom¬ 
etry measured from NIRC 2 observations. Best-fit prim ary (red 
line; LHS 2924, data from IBurgasser fc McElwainll200^ and sec- 
onda ry templates (blue h i^^i 2MASS J04070885-I-1514565, data 
from IBurgasser et al.l20Q4l) are shown at their relative scaling. The 
gray bars at top indicate the regions over which the fitting was 
done. 


scribed in B15. The NIRSPEC spectrum (Figure [7l) is of 
high quality (median S/N = 82) and similar to data re¬ 
ported in B15. The Hamilton data had a somewhat lower 
signal-to-noise (S/N) than prior observations (S/N k, 5 
at 7500 A) due to a shorter total integration of 3000 s. 

Both spectra were analyzed as described in B15 for 
radial motion and, for the Hamilton data, Ha emis¬ 
sion equivalent width (EW). Despite the low S/N, cross¬ 
correlation of the Hamilton spectrum with contempora¬ 
neous observations of the RV standard GJ 251 yielded 
an RV = -1-86.4±0.5 km s“^, significantly different than 
the mean motion reported in B15 (-|-82.5±0.4 km s“^). 
Marginal Ha emission was observed in these data, with 
EW = — 2±1 A, consistent with the lowest emission 
states observed in B15. 

For the NIRSPEC data, we re-analyzed the obser¬ 
vations reported here and the 2014 observations us¬ 
ing an updated Markov Chain Monte Carlo (MCMC) 
adaptation of the f orward-modeling method described in 
iBlake et’aTI (1201011 and B15. We us ed the BT-Settl at¬ 
mosphere models (lAllard et al.ll2012ll with u pdated solar 
abundance values from ICaffau et al.l (|2011[1 over an ef¬ 
fective temperature (Teff) range of 1600-2900 K and 
surface grav ity (logq) range of 4.5-5.5 ( cgs). The So¬ 
lar atlas of iLivingston fc Wallace! (| 199 111 was used to 
model telluric absorption features. Table [3] summarizes 
the RV and rotational velocities (usini) inferred from 
these analyses, while Figure 0 shows the best-fit model, 
with Teff = 2700 K and logg = 5.5, to the most recent 
NIRSPEC observations. The rotational velocity mea¬ 
surements are mutually consistent, with a mean value 
of ?;sinf = 8.6±0.4 km s“^. However, like the Hamil¬ 
ton observations, the most recent NIRSPEC RV mea¬ 
surement of -1-85.2±0.2 km s“^ is significantly distinct 
from the mean RV reported in B15 (-|-83.7±0.4 km s“^) 
and follows the trend of increasing recessional motion. 
With these independent measurements, we conclude that 
WISE J0720—0846A is an RV variable gravitationally 
perturbed by its brown dwarf secondary. 


Table 3 

Radial and Rotational Velocities from NIRSPEC Observations 


UT Date 

MJD 

RV 

(km s“^) 

usini 
(km s“^) 

2014 Jan 19 

56676.00968 

-|-84.0±0.3 

8.8±0.7 

2014 Jan 20 

56677.00048 

-l-83.2±0.2 

9.3±0.8 

2014 Mar 10 

56725.71832 

+82.9±0.5 

6.6±2.3 

2014 Apr 12 

56758.74056 

-|-84.3±0.4 

9.9±1.8 

2014 Dec 8 

56999.00802 

-l-85.2±0.3 

7.6±0.8 

Mean 


-l-83.8±0.8“ 

8.6T0.8 


^ = 35.9 indicates data inconsistent with a constant radial 


velocity. 


We note that the best-fitting models for the forward¬ 
modeling analyses were consistently in the ranges Tef / = 
2500-2800 K and logg = 5.0-5.5. The temperatures are 
somewhat higher than those expected f or an M9.5 dwarf 
(Teff ~ 2300 K; iSteohens et al.ll200^ . but the surface 
gravity is consistent with the lack of Li I absorption and 
low-surface gravity spectral features, indicating an age 
>1 Gyr. We verified that lower temperature and lower 
surface gravity models yielded identical RVs, so this mea¬ 
surement appears to be insensitive to the specific model 
over the parameter range examined here. A robust in¬ 
vestigation of the atmospheric parameters of this source 
is deferred to a later study. 


4. IMPROVED CONSTRAINTS ON THE ORBIT OF 
WISE J0720-0846AB 

4.1. Methodology 

With WISE J0720—0846AB verified as a gravitation¬ 
ally bound binary with orbital motion detected in all 
three spatial dimensions, we can begin to constrain the 
orbital properties of the system and the physical proper¬ 
ties of its co mponents. We adapted the MCMC analysis 
described in IBurgasser et al.l (1201211 to include both RV 
and relative astrometry measurements. We employed an 
orbit model with nine parameters, 

0 = (P,a,e,i,uj,U,Mo,q,VcoM,d) (1) 

where P is the period of the orbit in years, a is the semi¬ 
major axis in AU, e is the eccentricity, i is the inclination, 
u! is the argument of periastron, Q is the longitude of 
nodes, Mp is the mean anomaly at epoch tq = 2014.0896 
('MJEP^ = 56675.982), q = M2/M1 is the system mass 
ratio, VcoM is the center of mass (systemic) radial veloc¬ 
ity in km s“^, and d is the distance in pc. The primary 
radial velocity as a function of time t, Vi{t), is 

Vi{t) = Ki [ecosw -I- cos (T(t) -I- w)] -I- Vcom ( 2 ) 


where 


27rasinf q 

PVT^l + q 


(3) 


and the true anomaly T{t) is related to the eccentric 
anomaly E(t) through 


tan 


m 

2 


1 -I- e 
1 — e 


■ tan 


m 

2 


Modified Julian Date = Julian Date - 2400000.5 


( 4 ) 
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Wavelength (|j,m) 

2.29614 2.29971 2.30320 2.30660 2.30992 2.31317 2.31636 2.31947 



Pixels 


Figure 7. Extracted high-resolution (A/AA = 20,000) K-band spectrum of WISE J0720—0846 obtained with NIRSPEC on UT 2014 
December 8 (black line), compared to a best-fit model combining a Teff = 2700 K, logg = 5.5 atmosphere model from I Allard et II2012I . 
red line) with scaled telluric absorption (green line). The difference between data and model (O-C) is shown in black at the bottom of the 
plot and is dominated by fringing residuals; the ±1(7 uncertainty spectrum is indicated in grey. 


which is solved by Kepler’s Equation; 

M{t) — Mq = = E(t) — e sin E(t). (5) 

The angular separation vector from primary component 
to secondary component, p = (AQ;(t), A5(t)) is deter¬ 
mined from 


Aa{t) = 
ASit) = 



A{cosE{t) 

B{cosE{t) 


— e) ± E\/l — sin E{t) 

— e) + G\/l — sin E{t) 


( 6 ) 

(7) 


where Act and A^ are the angular separations on sky 
measured in arcsecond s, and A, B, F and G are t he 
Thiele-Innes constants (jlnneslIlQT^ Ivan den Boslll927[ ): 


= cos w cos n — sin 01 sin n cos J (8) 

B = cos 01 sin fl ± sin a; cos Id cos t (9) 

E = — sin o; cos Sd — cos 01 sin n cos z (10) 

G = — sinwsinfd ± cosoi cosfdcosz. (11) 


Note that the total system mass {Mtot = a^/P^ in solar 
masses) and component masses (Mi = Mto(/[l—g], M2 = 
gMi) could in principle be uniquely inferred from these 
parameters if sufficiently constrained. 

We selected an initial parameter set that visually co¬ 
incided with the observations through manual exper¬ 


imentation. We then computed a chain of 10^ pa- 
rameter sets using the Metropol i s-Has tings algorithm 
(|MetroDolis et al.l 1195.31 : iHastingsl 119701 ). at each step 
varying parameter 6j —f 0' by drawing a random offset 
from a normal distribution 


P{9j\9j)(xe 


( 12 ) 


where /3 is the set of jump stepsEB We applied addi¬ 
tional parameter constraints of 0.5 yr < P < 30 yr, 
e < 0.8 and 4 pc < d < 8 pc to eliminate improba¬ 
ble regions of parameter space; note that the eccentric¬ 
ity cutoff is be yond the e ~ 0 . 6 lim it suggested in em¬ 
pirical data by iDumiv fc Lid (|2011[ ). We also limited 
the component masses to 0.055 Mq < Mi < 0.15 Mq 
given the spectral classification of the primary and lack 
of Li I absorption in its optical spectrum (B15; 115), and 
M2 < 0.075 M0 given the substellar nature of the sec¬ 
ondary. Orbit models were compared to the data using 
a statistic that combined both RV and relative astro- 


11 We used in initial set /3 = (3 yr, 0.5 AU, 0.2,10°,10°,10°,10°, 
0.2, 1.0 km s“^, 1.0 pc). 
















































metric measurements: 


2 


x-E 
2=1 
N^st 

E 

i=i 






(model) \ 2 


Aaj' 


g (A^ 

i-1 


(o&s) 




(model) \ 2 




(13) 


where Nf/y and Nast are the number of RV and astromet¬ 
ric measurements, respectively, and cr the measurement 
errors. Model values were calculated at the same epochs 
as the observations. The criterion to adopt successive pa¬ 
rameter sets was t/(0,1) < where C/(0,1) 

is a random number drawn from a uniform distribution 
between 0 and 1. We compared models to data after each 
individual parameter change rather than changing all pa¬ 
rameters, a procedure we found greatly improved the ac¬ 
ceptance rate, which declined from 10% to 1% through 
the chain. Separate fits were made to the NIRSPEC -I- 
N1RC2 and Hamilton -I-N1RC2 measurements, and to all 
of the data. Note that in the last case we did not take 
into acco unt a poss ible velocity offset between the RV 
datasets (jFordll2nn^ . 

To test for convergence, we monitored parameter au¬ 
tocorrelation and the evolution of subchain variance 
through the chain. We also generated M = 10 chains of 
length 2n = 10® steps on the Hamilton -|- N1RC2 data, 
varying the initial parameter set Omit as 

O^mt =9p + N{0, l)ap (14) 

where 6p and (jp are the median and half quantile 
ranges of the posterior parameter distributions (see be¬ 
low) and N{0, 1) is a normal distribution centered on 
zero with unit variance. We quantified the within chain 
and between chain varia nce of all parameters using the 
iGelman fc RubinI (jl992tl scale reduction factor 

where Wk = Cfe is the within chain variance for pa¬ 
rameter Ok, equal to the average of parameter variances 
^jk = ;r^Z)”=i(%fc - Ojk)^ for each of the ten chains 
using the second half of the chains as the sample; and 
Bk = ~ ^kY is the between chain vari¬ 

ance, equal to the variance in parameter averages Ojk- 
Scale reduction factors were within 5% of unity for all 
modeled parameters with the exception of P (1.35), a 
(1.37) and uj (1.51). As described below, these three pa¬ 
rameters were weakly constrained, and this analysis sug¬ 
gests that the primary MCMC chains may not converge 
for our limited datasets, re-emphasizing that the orbital 
results presented here should be considered preliminary. 


4.2. Results 

Figure |S] shows the best-fit orbits from our separate 
analyses of the RV and imaging datasets. Table |4] lists 
the best-fit orbital and component parameters, as well as 


the median and 16% and 84% quantiles of the parame¬ 
ter distributions, after eliminating the first 10% of each 
chain. Figures |9] through [TT] display the distributions and 
correlations of parameters P, a, e, i, q and yitot for all 
three datasets. 

The best-fit orbits are exceptionally good fits for the 
individual RV datasets, with = 5.83 for nine data 
points (zero degrees of freedom) for NIRSPFC -I- NIRC2 
and = 2.84 for twelve data points (three degrees of 
freedom) for Hamilton -|- N1RC2. Note that these low y^ 
values largely reflect the underconstained nature of the 
solution in orbital phase; i.e., a relatively large family 
of solutions is formally consistent with these data. The 
combined dataset yields a somewhat poorer best ht, with 
modest disagreement between contemporary NIRSPEC 
and Hamilton RVs driving the y^ values. Nevertheless, 
this solution is nearly identical to that of the Hamilton 
-I- NIRC2 dataset. While the best-fit solutions between 
the NIRSPEC -h NIRC2 and the Hamilton -h NIRC2 
datasets are distinct, the median parameters are con¬ 
sistent within the uncertainties for all datasets. Since 
parameters inferred from the Hamilton -|- NIRC2 data 
are best constrained, we refer to these in the following 
discussion. 

Several of the parameters are very well constrained, 
most notably the orbital inclination (i = 93?6'®^ip, the 
longitude of ascending node (H = 82?7±2?1) and the cen¬ 
ter of mass radial velocity {Vcom = 4-87.4±0.8 km s“^). 
The tight constraints on i and H stem from the near- 
radial motion of the secondary between the two NIRC2 
epochs. The position angle of the secondary changed 
by only 5° ±2° between these two observations but the 
secondary has moved 40% further away; this is possible 
only when the orbit is viewed very close to edge-on or 
e ~ 1. As discussed below, even with an eccentricity 
at the proscribed limits, the inclination remains close to 
90°. Similarly, H is constrained by the orientation of the 
orbit on the sky. The tight constraint on Vcom arises 
from the small variance in observed RVs and consistency 
in the values measured in the earlier epochs. The NIR¬ 
SPEC -I- NIRC2 and Hamilton -|- NIRC2 datasets yield 
nearly identical values for all three of these parameters. 

Parameter distributions for P, e and q are more weakly 
constrained, and the latter two abut the imposed limits, 
emphasizing that we do not yet have sufficient coverage of 
the orbit of WISE J0720—0846AB to robustly determine 
them. To gain some insight on how the quality of the 
fits vary as e changes, we performed additional MCMC 
chains on Hamilton -|- NIRC2 data using models with 
hxed e = 0.2, 0.4, 0.6 and 0.8. Figure and Table [S] 
display the results of this experiment. We found that 
all parameters remain consistent with the unconstrained 
MCMC analysis, although q and M2 values are modestly 
smaller for the e = 0.2. Parameter uncertainties are con¬ 
siderably larger for the e = 0.8 case, in part reflecting the 
lower values of viable solutions. It is clear that the e 
= 0.2 and e = 0.4 models are poor representations of the 
RV data, with best-fit values that are significantly 
worse and can be eliminated with >95% confidence as 
compared to the unconstrained modelQ This analysis 

This confidence level was computed using the F-test proba- 
bility distribution function, comparing the unconstrained best-fit 
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Date (year) RA Separation (mas) 


Figure 8. Best-fit orbit from MCMC analysis based on NIRSPEC -|- NIRC2 data (top), Hamilton Spectrograph -h NIRC2 data (middle), 
and all data (bottom). The left panels show the radial motion of both the primary (solid black line) and secondary (red dashed line) motion 
compared to primary RV measurements (open circles with error bars). The right panels show the orbital motion of the secondary (blue 
line) relative to the primary (black dot at the origin) projected on the sky, compared to NIRC2 measurements (open circles; error bars are 
small on this scale). The arrow indicates the direction of orbital motion at apoapse. Note the different scales for the different fits. 


supports an eccentic orbit for WISE J0720—0846AB, but 
a firm constraint on its value should emerge with further 
monitoring. 

5. DISCUSSION 

Our follow-up observations of WISE J0720—0846AB 
verify both the magnetic and binary nature of this nearby 
system, and provide first constraints on the physical 
properties of the components. While these first mea¬ 
surements of orbital motion do not robustly constrain 
system or component masses, they do constrain the or¬ 
bital geometry, most notably inclination, which is within 
a few degrees of edge-on. The scale of this system is 
such that it is unlikely to eclipse; nevertheless, the or¬ 
bit orientation allows us to make some constraints on 
the rotational properties of the primary. The orbital and 

model (x^ = 2.84, 3 degrees of freedom) to the constrained models 
(x^ = 39.3 and 25.1, 4 degrees of freedom). 


primary rotational axes of main sequence binaries with 
close separations (<20AU) are generally aligned (|Weisl 
11974 iHald 119941) . The only very low mass binary for 
which spin-orbit alignment has been tested, 2MASSW 
J0746425-I-200032AB, shows alignment between orbital 
and rotational axes of both components to within 5° 
(jHardinsf et al.l[^013fl . Assuming similar alignment of the 
rotational and orbital axes of WISE J0720—0846A and 
a radius of 0.1 Rq, our vsini measurement implies a ro¬ 
tation period of 14.2±0.7 hr. This value is remarkably 
close to the marginally-indicated variability period from 
white light monitoring reported in B15, 14.00±0.05 hr for 
a 1.3±0.5% variability amplitude. The rotation period 
is somewhat s lower than the m ean for comparably clas¬ 
sified sources (|Irwin et al.ll20UT) and cons iderably slower 
than most highly variable L and T dwarfs (jRadigan et al.l 
I20l4 iMetchev et al.l I2015D . Its slow rotation rate and 
orientation may explain the weaker than average mag- 
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Table 4 

Parameters from Orbital Analysis 



NIRSPEC-I-NIRC2 

Hamilton-|-NIRC2 


All Data 

Parameter 

Best-tit 

Median 

Best-tit 

Median 

Best-tit 

Median 


Modeled Parameters 

Best 

5.83 


2.84 


27.6 


p- (yr) 
a (AU) 

i n 
c n 
n n 

Mo n 

g 

VcOM (km s“l) 
(pc) 

10.4 

2.4 

0.79 

92.1 

109 

81.7 

5 

0.36 

87.2 

4.5 

fi 1 

1 y-|-0.9 
-*-•'-0.5 

n yQ + 0.06 

'^-0.09 

oq-|-2d 

od_24 

82 . 2 / 2 ? 

20 /?? 

o- 6 i/g:?? 

87.3/°;® 

6.7/°;« 

2.9 

1.1 

0.79 

92.8 

70 

82.4 

10 

0.61 

87.3 

4.5 

4 1+2-'^ 

+•1-1.3 

1 q-tO.5 

0 77+0-02 
'-0.04 

93 . 6 /?? 

76/?? 

82.7//? 

13/1 

0 . 77 //?? 

87.4//j 

5-4/g;? 

3.1 

1.1 

0.80 

94.9 

54 

82.2 

24 

0.91 

87.4 

6.8 

4 5+2.7 

+•^-1.4 

1 4+0-5 
-*-•+-0.3 

0.76/°;°? 

94.4/?;« 

74-I-I6 
'+-16 
oq q”l"2.2 

oo.o_2 2 
17/? 

0.75/°;?° 

87.5/°;? 

6-2//,? 

interred Parameters 

Mtot (^Iq) 

Ml (Mq) 

M2 (Mq) 

Ki (km s“^) 

K 2 (km s“l) 

0.13 

0.098 

0.035 

3.0 

8.4 

0.080/°;°?i 
n o'si 

U.UOi_Q Q24 

4 4+1-2 
y q-I-2.2 
+ '^-1.6 

0.17 

0.105 

0.064 

7.2 

11.9 

0-15//^? 

0.082/°°?« 

0.064/??;^ 

6-3/?;J 

8-5/?-i 

0.16 

0.081 

0.074 

8.7 

9.5 

0-15//^? 

0.081//°?® 

0.062//°°? 

6-o//i 

8-2/?-i 

® Parameter was 

constrained to 

a limited value range 

in MCMC analysis. 







Table 5 




Parameters from Orbital Analysis 

of Hamilton -|- NIRC2 Data with Eccentricity Fixed 


Parameter 

e = 0.2 

e = 0.4 

e = 0.6 

e = 0.8 

Best X® 

39.3 

25.1 

10.4 

2.8 

(yr) 

4 4+0-8 

+•+-0.7 

4 3+1-0 
+•'^-0.8 

+•^-0.9 

3.8/?;? 

a (AU) 

1.40//?? 

-I qQ-|-0.23 

l.dO_Q 

1 4+0-3 
-*-•+-0.2 

1 q-|-0.6 

^•'^-0.3 

i n 

94.1/?;? 

94.4/?;? 

94.2/?;? 

93.5/?? 


0.69///® 

0.80//“ 

0.82//?? 

0.73//?° 

VcOM (km s-1) 

85.9//? 

86.5//? 

87.0//j 

87.3//,® 

^tot {^0) 

0-1+/.TO 

0.14//??? 

Q -I 5-1-0.02 
^•-*-^-0.02 

0-l+ao3 

Ml (Mq) 

0.083//°?? 

0 07Q+0-024 

u.uiy_0 014 

u.u/y_o 013 

0.083//??? 

M2 (Mq) 

0.059//???? 

'^•^'^+-0.011 

U.UOO_Q QJQ 

0.062//°?? 


^ Parameter was constrained to a limited value range in MCMC analysis. 


netic emission of WISE J0720—0846A, in both optical 
and radio bands, and its somewhat low flaring rate com¬ 
pared to other late M dwarfs. In addition, with a 1-2 AU 
semimajor axis, magnetospheric interaction between pri¬ 
mary and secondary is unlikely to play a role in d riving 
magnetic emission in either source (|Schriivei]l2009ll . 

A more promising role for WISE J0720—0846AB is as a 
testbed for brown dwarf evolution ary models. As demon¬ 
strated by several studies (e.g., ILiu et all [Mil IMnI: 
IKonopackv et al.l 120101 : IDuduv et al.ll201 Jb the inferred 
masses and atmospheric parameters of several very low 
mass binaries diverge to varying degrees from evolution¬ 
ary model predictions. Limiting factors for such analyses 
include lack of knowlege of detailed spectra for individual 
components, inability to measure component masses, low 
quality distance determinations, and/or lack of indepen¬ 
dent age determinations. While the last factor may be 
challenging to overcome for WISE J0720—0846AB, the 
proximity and detection of reflex motion in both astro- 
metric and radial coordinates allows the first three to be 
addressed, and we anticipate that this system will pro¬ 


vide a high-quality test of models in the next 5-10 years. 

Finally, we note that despite being able to detect 
reflex motion in the primary, the inferred close pas- 
sage of WISE J 0720— 0846AB to the Sun deduced by 
iMamaiek et al.l (j2015| l on the basis of the system kine¬ 
matics is not ruled out. The systemic motion estimated 
here is about 4 km s“^ larger than that assumed by 
these authors, so the impact parameter is roughly 5% 
smaller. The astrometric perturbation of the primary 
induced by the secondary, of order 100 mas, is suffi¬ 
ciently slow to play little role in modulating the parallax 
or proper motion of the system significantly given cur¬ 
rent measurement uncertainties. Improved monitoring 
of the systemic orbital motion will be needed to make a 
more precise estimate of the geometry and timescale of 
WISE J0720—0846AB’s closest approach, and the cor¬ 
responding perturbation it has made on our Oort cloud 
comet population. 

The authors thank Randy Campbell, Heather Her- 
shley, and Marc Kassis at Keck Observatory; and 
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Figure 9. Parameter distributions and correlations (triangle plot) 
for period (P), semi-major axis (a), eccentricity (e), inclination (i), 
mass ratio (q) and total system mass (Mtot) based on our MCMC 
orbital analysis for the NIRSPEC -h NIRC2 data. The fits assume 
a weak constraint on eccentricity (0 < e < 0.8) and period (0.5 yr 
< P < 30 yr). Contour plots show two-dimensional distribu¬ 
tions for parameter pairs, highlighting correlations. Normalized 
histograms at the ends of rows are marginalized over all other pa¬ 
rameters. Median values are indicated by solid lines in all panels; 
16% and 84% quantiles are indicated by dashed lines in the his¬ 
tograms. Imposed parameter limits for P and e are indicated by 
dotted lines. 


5 



Figure 10. Same as Figure 0 for the Hamilton Spectrograph -h 
NIRC2 data. 
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Figure 11. Same as Figure 0 for all RV and imaging data. 
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